Spin flip lifetimes in superconducting atom chips: BCS versus Eliashberg theory 
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We investigate theoretically the magnetic spin-flip transitions of neutral atoms trapped near a 
superconducting slab. Our calculations are based on a quantum-theoretical treatment of electro- 
magnetic radiation near dielectric and metallic bodies. Specific results are given for rubidium atoms 
near a niobium superconductor. At the low frequencies typical of the atomic transitions, we find 
that BCS theory greatly overestimates coherence effects, which are much less pronounced when 
quasiparticle lifetime effects are included through Eliashberg theory. At 4.2 K, the typical atomic 
spin lifetime is found to be larger than a thousand seconds, even for atom-superconductor distances 
of one micrometer. This constitutes a large enhancement in comparison with normal metals. 
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I. INTRODUCTION 

Over the last few years, enormous progress has 
been made in magnetic trapping of ultracold neutral 
atoms near microstructured solid-state surface, some- 
times known as atom chips [H, [3, H, 0] ■ The atoms can be 
manipulated through variation of the magnetic confine- 
ment potential, either by changing currents through gate 
wires mounted on the chip or by modifying the strength 
of additional radio-frequency control fields. These exter- 
nal, time-dependent parameters thus provide a versatile 
method of atom manipulation, and make atom chips at- 
tractive for various a ppli cations, including atom interfer- 
ometry 

d 1 0, i Wm . quantum gates [H [11 [H, [11 

and coherent atom transport [l5| . In addition, the atoms 
may be used as a sensitive probe of the electromagnetic 
properties of the surface in the neV (MHz) energy range. 
For ex amp le they can image surface currents in a normal 
metal |16j or vortices and flux noise in a type-II super- 
conductor [l3|. 

On the other hand, the proximity of the ultracold 
atoms to the solid-state structure, introduces additional 
decoherence channels, which limit the performance of the 
atoms. Most importantly, Johnson-Nyquist noise cur- 
rents in the dielectric or metallic surface arrangements 
produce magnetic-field fiuctuations at the positions of 
the atoms. Upon undergoing spin-flip transitions, the 
atoms become more weakly trapped or are even lost from 
the microtrap (isl . [Toj . Typically, the spin-flip transition 
frequencies for magnetically trapped alkali atoms are in 
the sub-MHz range, and the radiation-atom coupling is 
therefore strongly enhanced by being in the near field 
regime [lO, HH, Ull . For atom-surface distances of the or- 
der of one micrometer, the atom lifetime typically drops 
below one second, which constitutes a serious limitation 
for atom chips. It was shown in Ref. (2^ that in order to 
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reduce the spin decoherence of atoms outside a metal in 
the normal state, one should avoid materials whose skin 
depth at the spin-fiip transition frequency is comparable 
with the atom-surface distance. For typical experimental 
designs using metals such as copper or gold, however, the 
atom-surface distances are precisely in this range [l^ [3 ■ 

Superconductors could reduce the magnetic noise level 
significantly and thereby boost the spin flip lifetimes 
by many orders of magnitude. Indeed, superconduct- 
ing atom chips have already been fabricated and tested 
(23I [2I with the aim of realizing controllable composite 
quantum systems. Previous estimates of the lifetime en- 
hancement relative to a normal metal surface have given 
factors of tens [l^ or millions [2^, depending on the 
theoretical approach. Scheel et al. considered the 
energy dissipation in the superconducting state result- 
ing from the modified quasiparticle dispersion, whereas 
Skagerstam et al. [251] considered the screening of the 
current fluctuations by the superconductor. The two ap- 
proaches are difficult to compare, since they ignore the 
strong modification of either the imaginary part [2^ or 
the real part [2^ of the optical conductivity in the su- 
perconducting state. The question of how to describe the 
problem properly led to some dispute [1^ [131 . 

In this paper, we resolve the dispute and present a 
scheme for the proper description of magnetic spin-flip 
rates in atoms on a superconducting atom chip. Our 
analysis is based on three descriptions of superconduc- 
tivity. We start with the two-fluid model [25|, and then 
progress via the Bardeen-Cooper-Schrieffer (BCS) the- 
ory |28| to a more elaborate framework, the Eliashberg 
theory [29| , which we find is needed for a proper descrip- 
tion of this problem. For typical spin flip frequencies on 
a chip (IkHz-lOMHz), we point out that the BCS the- 
ory significantly overestimates the optical conductivity 
and hence gives too high a value for the spin flip rate. 
A realistic calculation of the conductivity [30*1 requires 
further elaboration, in the framework of the Eliashberg 
theory, to include lifetime effects of the quasiparticles due 
to phonon scattering. This results in a reliable estimate 
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of the spin flip rate, which ends up not far from the two- 
fluid result of Skagerstam et al. |;25j]- We conclude that 
superconducting surfaces can be used to achieve low spin 
flip rates in an atom chip, with lifetimes exceeding a thou- 
sand seconds for Rb atoms at 1 ^m from a Nb surface at 
4.2 K. 

We have organized our paper as follows. In Sec. II we 
present the three models for the description of supercon- 
ductors. Although there is already a vast literature on 
superconductivity, including many textbooks [Sll. [s^. [s^ . 
we give next a brief account of these approaches, mainly 
to make the paper as self-contained as possible. We dis- 
cuss the basic assumptions of the two-fluid model, the 
appearance and shortcoming of the coherence peak in 
BCS theory, and the description of quasiparticle damping 
and formation within the framework of Eliashberg the- 
ory. In Sec. Ill we present our results for the lifetime of 
an atom placed in the vicinity of a scmi-infinitc niobium 
sample. We compare the different approaches and discuss 
their respective advantages and disadvantages. Finally, 
we summarize our results in Sec. IV. 



II. THEORY OF SUPERCONDUCTIVITY 
A. Two-fluid model 

The first successful attempt to account for the elec- 
tromagnetic properties of superconductors was due to F. 
and H. London They devised a phenomenological 

two-fluid model that was able to explain many of the 
phenomena observed in superconductors. 

Within this model one assumes that there are two 
types of charge carrier, superconducting and normal, 
which react dilferently to external electromagnetic fields. 
We write n„(T) and ns{T) to denote the electron num- 
ber densities in the normal and superconducting states 
at temperature T, with n„(T) -I- ns{T) = no assumed to 
be constant. Although it does not become obvious from 
the two-fluid model itself, the superconducting carriers 
have to be associated with Cooper pairs. At temper- 
atures above the superconductor transition temperature 
Tc, only normal carriers are present and n„{T>Tc) = no, 
while at zero temperature all carriers are in the super- 
conducting state, ns{0) — hq. 

For the normal electrons, the response to a sufRciently 
weak external electric field E is given by Ohm's law 
jri = CriE, with in being the current density of the nor- 
mal electrons and (T„ the normal-state conductivity. For 
the superconducting current js, the London brothers in- 
troduced a new relation 



A 



dt 



E. 



(1) 



where A is a constant whose value varies for different 
superconducting materials. This describes the dynam- 
ics of carriers that are accelerated freely in an electric 
field. For a superconductor made up of free electrons (or 



indeed of free Cooper pairs), the value of A would be 
m/ (rise^), where m and e are the single electron mass 
and charge. In fact this also provides a useful estimate 
for real superconductors. Later in the paper we will re- 
write this relation in terms of the plasma frequency Wp, 
as A ~ l/(eoWp). As a consequence of the London equa- 
tion, ID), a static magnetic field can only penetrate into 
a superconductor by a distance of order Al — {A/ fio)2 
[sH . For this reason is called the penetration depth 
or London length. 

Consider an electric field oscillating as exp(—iujt). The 
response of the superconductor is given by 



jA 



E. 



Here, the expression in parentheses 



(2) 



(3) 



is known as the optical conductivity, though in this paper 
we will be using it at radio frequencies. We have intro- 
duced the skin depth S = (2/ fiQuan)'^ associated with 
the normal charge density. 

For the two-fluid model, Eq. ^ can be further sim- 
plified by noting that the two contributions vary with 
temperature only through the normal and superconduct- 
ing charge densities. Thus, with ctq being the conductiv- 
ity in the normal state and Aq the A-parameter at zero 
temperature, we have 



nn (T) 



no 



ujAo 



1 



no 



(4) 



For T < Tc, a suitable form for the temperature depen- 
dence of the normal density is provide d by the Gorter- 
Casimir expression nn{T) — {T/TcYno [35|. 



B. Bardeen— Cooper— Schriefffer (BCS) theory 

Despite its success, the London theory has a number 
of shortcomings. First, it is phenomenological and not 
based on a microscopic model. Second, its predictions 
cannot account for all experimental observations. A rel- 
evant example here is its inability to account for the so- 
called coherence peafc, that was first observed in NMR by 
Hebel and Slichter [3g| . This peak is most pronounced at 
low frequencies and is thus of importance for the analysis 
of spin decoherence in superconducting atom chips. In 
order to understand its origin we introduce the theory of 
Bardeen, Cooper, and Schrieffer (BCS) [H. 



1. BCS ground state 

BCS theory is based on Frohlich's observation [s^l 
that electrons close to the Fermi energy tp can attract 
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each other through the exchange of virtual phonons, and 
Cooper's demonstration ^3^] that due to this interaction 
the Fermi sea is unstable against the formation of a cer- 
tain kind of quasi-bound pair. The attractive electron- 
electron interaction is usually described by the pairing 
Hamiltonian [3l| 

kcr k,k' 

Here cj^.^ is the field operator for the creation of an elec- 
tron with wavevector k and spin orientation cr, = 
Ek — ef is the single-electron energy eko- measured with 
respect to ep, and V is the strength of the attrac- 
tive phonon-mediated electron-electron interaction. The 
prime on the sum indicates that this interaction has to 
be considered only for electrons with energy smaller than 
the Debye energy Hujd- 

As result of this coupling, electrons are promoted from 
states below the Fermi energy to states above to form 
Cooper pairs. This process comes to a halt when the in- 
crease in kinetic energy is no longer compensated by the 
reduction in potential energy from the pairing. To model 
the phase transition associated with the formation of 
Cooper pairs, one assumes that the interaction operator 
c_^j^Cj.| is practically a c-number b'^, with small fluctua- 
tions about this value. One then formally writes all pairs 

of operators in the form c_j,|C^| = ^k + ('^-fci'-kT ~ ^k) 
and neglects the terms bilinear in the parenthetical quan- 
tities. The resulting mean-field Hamiltonian can be di- 
agonalized through a Bogoliubov transformation 

CkT = "k7ko + '"k7ki ' cLkx = ~"k7ko + "k7ki > 

where 7^q and 7^^^ create Fermionic quasiparticles that 
are linear superpositions of the bare electron states, and 
the coefficients Uk and fk are chosen to diagonalize the 
Hamiltonian, 

Hbcs = XI (7ko7ko + 7ki7ki) + const . (6) 

k 

Here, = (Ck + A^)^ are the new quasiparticle exci- 
tation energies in the superconducting state, and A — 
V'^^'b'^ is the order parameter or gap parameter. A 
has to be determined from the numbers 6k which are the 
thermal and quantum averages of 

bl = tr (e-^^--c_k^Ck^) /tre-'^^Bcs ^ (7) 

where (3 = l/{kBT). Equation ([7]) is a self-consistency 
relation, since the values of 6^ are hidden within i?Bcs 
through its dependence on the quasiparticle energies E^. 
Thermally excited quasiparticles with energy i?k restrict 
the phase space available for forming Cooper pairs and 
thereby reduce the gap parameter A. 



2. Coherence peak 

The density of these quasiparticle states at energy E 
is given by [31| 

[0, < A. 

At zero temperature, no quasiparticles are excited and 
therefore the only way to deposit energy in the super- 
conductor is to break up Cooper pairs. Consequently the 
real part a' of the T = conductivity is strictly zero for 
electric field frequencies below 2A/h. At non-zero tem- 
peratures however, many quasiparticles may be excited 
just above the gap because the density of states is so 
high there — indeed p{E) diverges in Eq. ^ at E = A. 
This opens up a mechanism for dissipation at low fre- 
quency. The corresponding a' involves the density of 
quasiparticles, which is proportional to p{E), and the 
density of final states for absorption of a photon at fre- 
quency uj, which is proportional to p{E -\- Huj). Integra- 
tion over E produces a logarithmically divergent conduc- 
tivity o''(a;) (To lii(2A/u;). This enhancement, which 
was first observed in nuclear magnetic resonance [s^, is 
known as the Hebel-Slichter or coherence peak. This rea- 
soning is supported by Mattis and Bardeen's expression 
for the optical conductivity [39], which was computed 
with the random-phase approximation and in the dirty 
limit, where scattering by impurities reduces the coher- 
ence length to less than the magnetic-field penetration 
length Xl. This gives the same logarithmic divergence 
of (t'(ijj) at low frequency [32, 40]. At zero frequency, we 
note that a' has another singularity of S type, associated 
with the dc response of the superfluid. 

For the sub-MHz spin-fiip transitions of magnetically 
trapped ultracold atoms, the BCS theory thus predicts 
a strong modification of the optical conductivity in com- 
parison to the frequency-independent value of Eq. (|4]) 
given by the two-fluid model: a' = ao{nn / tiq) ■ 

C. Eliashberg theory 

While the BCS theory incorporates the mixing of free 
electron states through their coupling to virtual phonons, 
it does not include the dissipative effects associated with 
the emission and absorption of real phonons. This broad- 
ens the quasiparticle states and softens the divergence of 
the conductivity at low frequency so that it is much less 
dramatic. 

Phonon scattering converts the electron wavevector k 
to wavevectors k' at a rate l/rk, given by Fermi's Golden 
Rule as 

1 27r ■r-^ I A |2 r/ j- A\ 

— = 1- fk.k' "^(^k - ek' - ) 

k' 

X [2ntY,{nw^) + 1] , (9) 
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where hiu^ is the energy of a phonon in mode A with 
wavevector q = k — k', n^h is the number of thermal 
phonons in the mode and j^/ is the off-diagonal matrix 
element of the the electron-phonon interaction Hamilto- 
nian. Since electron energies are typically two orders of 
magnitude larger than the Debye energy, the phonon en- 
ergies entering the Dirac delta function in Eq. ^ can be 
safely neglected. This approximation leads one to define 
the dimensionless quantity 

a^F^iuj) = kk.k' r '5(ek' - eF)6{cu - u;^) , (10) 

k',A 

known in literature as the Eliashberg function ^33]. Thus, 
the electron scattering rate at low temperatures can be 
conveniently written as 

27r 

I/tw^— div a'^ Fkiui) [2n{ha;) + 1] . (11) 
" Jo 

In this expression, the Eliashberg function encapsulates 
all the relevant information about the electron-phonon 
coupling and the Fermi surface. The complex self en- 
ergy E resulting from this coupling gives both the scat- 
tering rate that we have just discussed, through fi./Tk — 
2 3m[E(ek)], and the energy shift 3fie[S(ek)] of the elec- 
tron. 

There are two kinds of self-energy function in the de- 
scription of a superconductor, usually labeled normal and 
anomalous. The normal component has the same mean- 
ing as in an ordinary metal, whereas the anomalous one is 
directly related to the opening of the gap due to the for- 
mation of Cooper pairs. These are closely related because 
the scattering rate and distortion of the electron bands 
due to the electron-phonon coupling depend strongly on 
the superconducting gap, and vice versa. This inter- 
dependence is accounted for by the so-called Eliashberg 
equations, which must be solved self-consistently [33| . 

A powerful numerical implementation for the solution 
of the Eliashberg equations has been dev elop ed by Car- 
botte, Marsiglio, and coworkers [30|, EH, lUSl], where 
one first computes the electron Green function in Mat- 
subara space and then performs an analytic continuation 
by means of an iterative procedure. The real-frequency- 
axis Green functions can be used finally to compute 
the optical conductivity [13, [H, El] , including not only 
electron-phonon interactions, described above, but also 
the effects of elastic impurity scattering. 

D. Impurity effects 

We conclude this section by briefly addressing effects 
due to elastic impurity scattering. In conventional su- 
perconductors impurities are deemed to be innocuous as 
a result of Anderson's argument [i^, El] , which goes as 
follows. In the normal state, the electrons can be de- 
scribed by wavefunctions (/)„|(r) and 0„j(r), where (/>„ 
is supposed to include the effects of impurity scattering. 



The quantum number n replaces the wavenumber k of 
the pure metal. In the pure superconductor, the Cooper 
pair is composed of the states (k, t) and (— k, J,). An- 
derson pointed out that the second of these states is the 
first with momentum and current reversed in time. In 
an impure superconductor the main contribution to the 
pairing should be also between the time-reversed states 
(/>„l(r) and (/)*|(r). The pairing Hamiltonian (O can thus 

be expressed in terms of the new operators and cl^^, 
where the interaction matrix element between two states 
becomes 

Vr.„,=Vj2\{n\k)\'\{n'\k')f = V. (12) 

k,k' 

Owing to the completeness relation of the states involved, 
the pairing Hamiltonian is not modified in the new basis 
(j)n- For this reason the superconductor properties such 
as, e.g., transition temperature, gap parameter, or quasi- 
particle density of states, are not significantly changed 
by the presence of impurities. 

The argument above applies not only to BCS but also 
to Eliashberg theory as long as the Eliashberg function 
Q!^Fk(w) has little dependence on the direction of k. 
This is indeed the case for the conventional s-wave su- 
perconductors we are considering. Moreover, any small 
anisotropy is randomized by the impurity scattering, so 
it suffices in this work to consider the average over all 
directions a'^F{u;) = (a^Fk(cLi)). 

Although impurities do not affect the pairing Hamilto- 
nian, the scattering from impurities at rate 7 plays an im- 
portant role in the electron transport because the normal 
conductivity ctq is approximately inversely proportional 
to 7. In the two-fluid model and in BCS theory, a'{uj) 
increases in direct proportion to (Tq as the scattering rate 
is reduced. In Eliashberg theory however, the situation is 
complicated by the presence of the inelastic phonon scat- 
tering, which tends to reduce the conductivity through 
the broadening of the density of quasiparticle states. As 
7 is reduced, this effect becomes relatively more impor- 
tant, causing (j'{lo) to increase more slowly than ctq- In 
the calculations that follow, we will allow 7 to be a vari- 
able in the optical conductivity [30, 43] so that we can 
explore this effect. We will find that this provides a con- 
nection between the two-fluid and BCS results as well as 
allowing us to make contact with real materials. 

III. RESULTS FOR THE ATOM TRAPPING 
LIFETIME 

We turn now to the spin flip rate for an atom located 
in vacuum near a superconducting slab, as illustrated in 
Fig. HI Following Refs. 21, 22, 25], we consider a ground- 
state alkali atom, magnetically trapped in a weak-field- 
seeking Zeeman sub-level. The noise in the magnetic 
field, due both to vacuum fluctuations and to thermal 
currents in the surface, induces transitions between the 
levels, making the atomic spin change direction (spin flip) 
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TABLE I: Normal-state conductivity ao measured on several 
different samples of niobium. The approximate scattering 
times r are obtained from the Drude model (|15p . The corre- 
sponding plasma frequency is ujp ~ lOeV/ft ~ 1.5 x lO^^s"^. 



Reference tJo {jiO. ^cm ^) r (fs) 

Perkowitz et al. [48] (X2 l0~~ 

Pronin et al. [49] 0.25 13 

Klein et al. [40] 0.85 43 

Casalbuoni et al. [47] 20 1000 



FIG. 1: (color online) Schematic geometrical setup. A plane 
metallic or superconducting slab lies parallel to the {x, y) 
plane. The atom with magnetic moment fi indicted by the 
arrow, is located in vacuum at a distance z from the surface. 
The atom suffers spontaneous or thermally stimulated mag- 
netic spin-flip transitions, as indicated by Go and G, thereby 
becoming more weakly trapped and eventually lost. Johnson 
current noise (jj) within the penetration depth A contributes 
to magnetic-field fluctuations at the position of the atom. 



and ultimately causing the atom to be lost from the mi- 
crotrap. 

As briefly outlined in Appendix \^ the spin-flip life- 
time of an atom at position ta is directly related to the 
imaginary part of the dyadic Green tensor G(rA, fyi, 'j-') 
of Maxwell's theory. The usual, free-space spontaneous 
emission rate is determined by the vacuum contribu- 
tion Go- For a typical transition frequency of /a — 
a;A/(27r) — 500 kHz, corresponding to an energy of ap- 
proximately 2 neV, this natural lifetime at zero temper- 
ature is To ~ 2 X 10^^ seconds [22l|, which can safely be 
considered infinite. The dominant contribution to the 
lifetime reduction comes from the magnetic-field fluctu- 
ations induced by the Johnson-Nyquist noise in the di- 
electric body. As shown in Fig. [1] and discussed in the 
Appendix, the current noise translates through the Green 
tensors to a magnetic-field fluctuation at the position of 
the atom. 

For a thick superconducting slab described by an op- 
tical conductivity in the limit (t"(u) ^ CT'(tj) and in the 
near field regime ^ z ^ 27r/fc one calculates, using 
the results of Ref. ^25], a spin- flip rate of 



— « Fo (nth + 1) 

TA 



1 



27 



64(tJ^o)l/2p^4 (cr")3/2 

(13) 

Here ta is the spin flip lifetime, Fq is the free-space 
decay rate, nth is the mean thermal photon number 
at the transition frequency to a, k — lua/c, and z is 
the atom-superconductor distance. In the following sec- 
tions we investigate the consequences for this rate of us- 
ing the expressions for the optical conductivity ct^lo) = 
a' {Lu)+ia" (uj) obtained from a two-fluid description, from 
BCS theory and from Eliashberg theory. 



A. Two-fluid model 

To estimate the order of magnitude of these parame- 
ters, we first consider the simple two-fluid model. Using 
the expression ^ for the optical conductivity, Eq. (fT3|) 
reduces to the expression [25| 



1 

TA 



To (nth + 1) 



1 



1 Ai 



(14) 



We take Nb as a representative superconducting material 
throughout. Table |T] shows a few values reported in the 
literature for the conductivity cto of the normal state. We 
note that the ultra-pure niobium sample of Ref. [i^l has 
a hundred times higher conductivity than the films of 
Refs. [ii,!!^. Through the simple Drude model ^ 



(15) 



we can relate do to an electron lifetime r = I/7 due to 
elastic scattering at impurities or defects, hujp is the bulk 
plasmon energy, which we set equal to 10 eV [H, [soj . 
The corresponding r values are given in the last column 
of Table HI With an atomic transition frequency of 500 
kHz, we obtain for the ultrapure sample a normal-state 
skin depth of So = y^2/(/ioW(7o) ~ 16 /im and a value 
approximately ten times larger for the films. 

A rough estimate for the penetration depth of the su- 
perconductor at zero temperature is given by 



A, 



1 



Mo £0 



c 



20 nm , 



(16) 



where we have assumed that all electrons move freely. 
This simple estimate is comparable to the BCS value of 
35 nm [51|, and to the experimental values of 46 nm for 
the ultrapure sample p7| and 90 nm for the niobium film 
in iia]- 



B. BCS versus Eliashberg theory 

Now we discuss how the two-fluid estimates are mod- 
ified within the framework of BCS and Eliashberg the- 
ories. For the BCS theory of niobium we use a zero- 
temperature gap parameter of A = 1.4 meV, correspond- 
ing to a transition temperature of Tc = 9.2 K, and a De- 
bye temperature of hLUu/kB = 275 K, and we compute 



6 



two fluid 

BCS 

A Y= 1 meV 
• 7= lOmeV 
Y= 100 meV 



e 1.5 




T/T 



FIG. 2: (color online) Temperature dependence of (a) real 
part a'{u!A) and (b) imaginary part a"{ujA) of the optical 
conductivity, normalised to the normal state conductivity (tq. 
ua = 2tt X 500kHz is the atomic spin-flip frequency, Tc — 9.2 
K is the superconductor transition temperature. The different 
lines correspond to the results for the two-fluid model (dashed 
line), using So — 16 /xm and Al(0) = 35 nm, BCS theory 
(solid line), and Eliashberg theory (symbols) for three elastic 
impurity scattering rates 7. 



the optical conductivity by means of the Mattis-Bardeen 
formulas in the dirty limit (soj . 

For the implementation of the Eliashberg equations, 
we have considered an a'^F{uj) function calculated using 
linear response theory [s^ and norm-conserving pseu- 
dopotentials. The electron-phonon matrix elements were 
calculated on a 32^ wavevector grid for both electrons 
and phonons. Our result (not shown) is similar to that 
presented in Ref. f53| , though with spectral features that 
are less pronounced, in better agreement with the data of 
tunneling experiments. As far as the calculated atomic 
spin flip rates are concerned, we do not find any signifi- 
cant difference between these two a^F(ijj) functions. 

Figure [5] shows results for the (a) real and (b) imagi- 
nary part of the optical conductivity versus temperature 
T. The solid lines show the results from BCS theory 
in the dirty limit, the dashed lines are for the two-fluid 
model and the symbol series are for Eliashberg theory 



with various values of the elastic impurity scattering rate 
7. In Fig. (Ha), the a'{coA) obtained from the two-fluid 
model decreases monotonically with temperature because 
of the decrease in the normal density ri„(T), whereas 
the BCS and Eliashberg curves show an enhancement of 
a'{LLjA) at temperatures immediately below the transi- 
tion temperature Tc. This is due to the coherence peak, 
which forms as a consequence of the modified quasiparti- 
cle dispersion in the superconducting state. The peak is 
most pronounced within the BCS framework in the dirty 
limit. As we move away from the dirty limit towards a 
clean superconductor, using the Eliashberg theory with 
decreasing rates 7, we observe that the peak gradually 
disappears, in agreement with [H^. Thus the Eliashberg 
theory interpolates between the two extreme cases of the 
two-fluid model and the dirty limit of BCS theory by 
varying the chosen value of 7. 

In Fig. ^h) we show the imaginary part of the optical 
conductivity. Again, the BCS result is according to the 
theory of Mattis and Bardeen 39] for a dirty supercon- 
ductor. Here too, we see that the Eliashberg theory with 
variable 7 provides a link between the two-fluid and BCS 
extremes. In the low frequency limit, the BCS result 
takes the analytical form 



,, , , ttA , a 

-BCsM-'^o^tanh^ 



(17) 



where A is the temperature-dependent gap parameter. 
We note that this 1/oj dependence of cr" is the same in 
all three models. This is a consequence of the Kramers- 
Kronig relations together with the fact that a' has a S 
singularity at = associated with the response of the 
superfluid to a dc field. 

We are using here a theory in which the material re- 
sponds locally to a field. Although this is not strictly 
so, nonlocality can be incorporated empirically into the 
theory of the superconductor through a modified pen- 
etration depth [2^. The effect of nonlocality on the 
atom-surface response is negligible since the penetration 
depth is small compared with the atom-superconductor 
distance. 

As discussed in Sec lIIB 2l the BCS coherence peak il- 
lustrated by the solid line in Fig. ^a.) increases with de- 
creasing frequency, diverging as w — > 0. This behaviour 
is greatly suppressed when inelastic phonon scattering is 
taken into account using Eliashberg theory (see also the 
discussion in Sec. IIIPp . as plotted in Fig. [3] This figure 
shows the real part of the optical conductivity (T'(a;) at 
4.2 K with three values of 7, spanning the 10 — 1000 ps 
range of scattering times given in Table [H The peak in 
Figure [3] at 1 THz is the conductivity associated with 
the breakup of Cooper pairs. The lower-frequency peak, 
which is the one of relevance here, no longer diverges at 
low frequency but reaches a constant value, shown in- 
set in the figure for frequencies below 2 GHz. Here, as 
in Fig. [21 the value of a' is normalised to the normal 
state conductivity to remove most of the dependence on 
7. Since the atomic spin flip frequency is bound to be 
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FIG. 3: (color online) o-'{(^) at 4.2 K as a function of frequency 
for three elastic scattering rates 7, as computed within the 
framework of Eliashberg theory. The peak at zero frequency 
is attributed to the condensate, and the peak at 1 THz to 
the breaking of Cooper pairs. In the inset we show that the 
condensate peak saturates at low frequencies. 



in this low frequency range, the Eliashberg results shown 
in Fig. m apply to all cases of experimental interest. We 
recall that the a' of the two fluid model in Eq. ([3]) is also 
frequency independent. 

Finally, in Fig. [3I^a) we show the spin-flip lifetime 
TA = 1/r [see Eq. ([T^] as a function of temperature 
for an atom-surface distance of 10 /im. The dashed line 
indicates the results obtained from the two-fluid model of 
Ref. [2^ . The lifetimes obtained from Eliashberg theory 
(symbols) are smaller, but only by a factor of ten or less: 
the influence of elastic scattering rates 7 on the spin-flip 
lifetime is not very strong. This indicates that the qual- 
ity of the niobium is not critical. Surprisingly, we find 
that ta is smallest for the high-quality sample with 7=1 
meV, highest for the intermediate value 7=10 meV, and 
falls off again slightly for 7 = 100 meV. In Fig. Hfb) we 
show TA as a function of atom-surface distance at 4.2 K 
(T/T, = 0.46). 

For an atom-surface distance of 1 /xm we obtain for 
7 = 1 meV a lifetime ta ~ 5000 seconds at a transi- 
tion frequency of (jj^/Stt = 500 kHz. Values for other 
distances can be obtained directly from the scaling of 
our central equation, Eq. (|13p. For other (low) frequen- 
cies, the lifetime given by Eq. p3p scales approximately 
as uja- This follows from the frequency independence of 
a' for uj > and the 1/uj dependence of a". 



IV. SUMMARY 

In this article, we have resolved the controversy sur- 
rounding the appropriate use of model assumptions for 
the electromagnetic energy dissipation in superconduct- 
ing materials. We have discussed the three most common 
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FIG. 4; (color online) Spin-flip lifetime ta of a trapped atom 
near a superconducting slab as a function of (a) temperature 
(for 10 /im atom-surface distance), and (b) distance (at 4.2 K). 
The smallest distance in (b) is 1 /^m (with ta ~ 5000 s for 
7 = 1 meV). For the calculation of ta we use Eq. (|13|l and 
the optical conductivity computed within Eliashberg theory 
for different elastic scattering rates 7. The atomic transition 
frequency is fixed at 500 kHz throughout. The dashed lines 
correspond to calculations performed with the two-fluid model 
and the parameters given in Ref. [2^ (London length Al = 35 



models of superconductivity, the two-fluid model, BCS 
and Eliashberg theories, in ascending order of sophisti- 
cation. The spin flip lifetime of neutral atoms trapped 
near a superconducting niobium surface is predicted to 
be much shorter when treated in the BCS theory than 
it is in the two fluid model. However, Eliashberg theory, 
which improves upon the BCS theory by including the 
finite quasiparticle lifetime, predicts only slightly shorter 
lifetimes. 

The Eliashberg theory interpolates between the two- 
fluid model and the BCS theory. For intermediate scat- 
tering rates, corresponding to real samples, the simple 
two-fluid model gives remarkably accurate estimates of 
the Eliashberg results. We have found that the lifetime 
depends only little on the precise value of the impurity 
scattering rate 7. 

Our numerical results based on the Eliashberg the- 
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ory show that the expected spin-flip Hfetime for an atom 
placed one micrometer away from a 4.2 K superconduct- 
ing planar niobium surface exceeds several thousand sec- 
onds at an atomic transition frequency of 500 kHz. This 
is expected to scale roughly as uj'\ and z^. Hence, su- 
perconducting surfaces provide an extremely low-noise 
environment for magnetically trapped neutral atoms and 
thus have great potential for coherent manipulation of 
atoms. 
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APPENDIX A: SPIN FLIP LIFETIME 

In this appendix we sketch briefly how to derive our 
basic expression ([13]). The derivation follows closely the 
general framework of Refs. [13, [Hj. There is, however, 
a subtle point regarding the fluctuation-dissipation the- 
orem, which we shall partly rephrase in the language 
of solid-state physics. In the general framework devel- 
oped by Welsch and coworkers [H, [H, one introduces 
Langevin noise operators with bosonic commutation re- 
lations, in order to fulfil the linear fluctuation-dissipation 
theorem. However, for calculating expectation values of 
bilinear operator products as in case of the spin-flip life- 
time, there is no particular need for such an approach. 

We consider an atom located in the vicinity of a di- 
electric body, as depicted in Fig. [U which is in a given 
magnetic sublevel. The coupling to the magnetic-field 
fiuctuations is described through a Zeeman interaction 
Hamiltonian in the rotating wave approximation. For a 
low-field seeking atom, the spin flip transition is associ- 
ated with an emission process, and the transition rate is 
simply given by Fermi's golden rule [1, 

a,/3 

(Al) 

Here a and /3 denote the Cartesian components, i and 
/ are the initial and final state of the scattering pro- 
cess, respectively, /i is the magnetic moment operator, 
and 'B{rA,oj) the Fourier transform of the magnetic field 



component with positive frequency. The position of the 
atom is ta and lua is the transition frequency. 

Now we relate the spectral density of the magnetic field 
to the current noise in the dielectric. In linear response 
theory we can use the Green tensor G of Maxwell's theory 
to relate the current j to the magnetic field B according 
to 

B(r,a;) = /.o J V x G{r,r' ,uj) j{r' ,uj) , (A2) 

with a corresponding equation for B'''(r,ti;). Thus, the 
spectral density of the magnetic-field fluctuations is given 
by convolving the spectral density (ja(r, w)j^(r', w)) of 
the current fluctuations with Maxwell's Green tensors, 
which describe how the held produced by the current 
fluctuation propagates to the position ta of the atom 
(see Fig. [ij . In the following we consider for simplicity 
only isotropic and local dielectric media. 

The calculation of (j(u))j'^(u!)) is a common problem 
in solid state physics [33|. For instance, in our present 
approach j could be the normal current j„ or the su- 
per current js of the superconductor. To express the 
spectral density of current correlations in terms of the 
optical conductivity, we first note that cr'(io) is related 
to the retarded current-current correlation via uj(7'{uj) — 
3fie /g°° <ii e*"* jUO)!}: which is a general result of 

linear-response theory [58l |. A common link between the 
ordered and retarded current-current correlation is pro- 
vided by the spectral function p{t) = (^[j (t) , {0)]) . Its 
Fourier transform can be obtained upon insertion of a 
complete set of states \m) with energy Em [1^ 

x2TTS{u;~{En~E„,)/h), (A3) 

with Z being the partition function. From Eq. (|A3p one 
immediately obtains {j{uj)j^uj)) = p(w)/(l — e^'"''^). 
The relation between the retarded current-current cor- 
relation and the spectral density is given through the 
Lehmann representation as p{uj) — 2uja'{uj) [33]. Thus, 
the desired relation between (j(cj)j'l'(a;)) and the optical 
conductivity reads 

(j(c^)jt(w)) = [nti^itiuj) + 1] 2uj<j'{ij) . (A4) 

One finally uses the expression 3m G = (po/uj) Ga'G* , 
which follows directly from Maxwell's equations fssl. [59|. 
to relate the scattering rate (lAip to the imaginary part 
of the Green tensor. The final equation is obtained 
according to the prescription given in Ref . [25j . 
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